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We study the closed Hamiltonian dynamics of a free particle moving on a ring, over one section of 
which it interacts linearly with a single harmonic oscillator. On the basis of numerical and analytical 
evidence, we conjecture that at small positive energies the phase space of our model is completely 
chaotic except for a single region of complete integrability with a smooth sharp boundary showing 
no KAM-type structures of any kind. This results in the cleanest mixed phase space structure 
possible, in which motions in the integrable region and in the chaotic region are clearly separated 
and independent of one another. For certain system parameters, this mixed phase space structure 
can be tuned to make either of the two components disappear, leaving a completely integrable or 
completely chaotic phase space. For other values of the system parameters, additional structures 
appear, such as KAM-like elliptic islands, and one parameter families of parabolic periodic orbits 
embedded in the chaotic sea. The latter are analogous to bouncing ball orbits seen in the stadium 
billiard. The analytical part of our study proceeds from a geometric description of the dynamics, 
and shows it to be equivalent to a linked twist map on the union of two intersecting disks. 



I. INTRODUCTION 



The one-dimensional free particle and the one-dimensional harmonic oscillator are arguably the two simplest quan- 
tum mechanical systems that exist. Nonetheless, and in spite of the intrinsic simplicity of such systems when treated 
in isolation, the problem of an essentially free particle interacting locally with one or more oscillators is difficult and 
important to a large class of physical systems. It arises in a variety of contexts, ranging from fundamental studies 
aimed at understanding the emergence of dissipation in Hamiltonian systems , to electron-phonon interactions in 
solids 0,13, and, more recently, to basic issues associated with the phenomena of quantum decoherence. In the con- 
densed matter literature, in particular, a great deal of theoretical work has focused on the nature of electron-phonon 
interactions, and the rich variety of behavior that occurs, such as the emergence of "polaronic" quasi-particles. In 
one version of the well-known Holstein Molecular Crystal Model, a tight-binding electron in a crystal moves between 
different unit cells, in each of which it interacts with a local oscillator Even in the simplest "spin-boson" form of 
this problem, in which the particle can be conceived as moving between just two sites, and in which it interacts with 
a single collective oscillator, the problem is not exactly soluble, and has been the subject of intense investigations 
regarding the appropriateness of various semi-classical approximations 4| . 

In this paper we show that the situation can be just as complex in completely classical versions of the problem. 
We study here the surprisingly rich classical dynamics of what is perhaps the simplest Hamiltonian model that one 
can think of that incorporates the essential features of this local free particle-oscillator interaction . Specifically, we 
consider a single classical particle of mass m, position x, and momentum po that moves on a ring, the circumference 
of which is divided into two sections. On one section the particle moves freely, on the other it interacts with a single 
oscillator of mass M and frequency oj. The associated Hamiltonian we write in the form 

i^^ll + ^ + iM-^X^-W.), (1) 

where P is the oscillator momentum, and Fq describes the strength of the interaction, which is linear in the oscillator 
coordinate X, but not in the particle position x. In the Hamiltonian (^), p(x) is a form factor localized around x — Q 
that describes the range of the interaction. We take p to equal unity throughout the interaction region |x| < ct, and 
to vanish on the section of the ring a -\- L/2 > \x\ > a lying outside this range. With this choice, except at those 
moments of the evolution when the particle arrives ai x = ±a, the particle and the oscillator are effectively uncoupled, 
and evolve independently. At "impacts", i.e., when the particle arrives at the edges of the interaction region, the 
particle receives an impulsive kick from the oscillator that conserves the total energy of the system. After each such 
kick, the particle again moves as a free particle either inside or outside of the interaction region. It is this essentially 
uncoupled evolution of the two subsystems between impulsive kicks of the particle at a; = ±(t that makes it possible 
to numerically and analytically track the dynamics. 

Despite its apparent simplicity, and as can be seen in the numerically determined phase plots appearing in Figs. 
I1I4I the model has a striking combination of features that are not often seen together in a single closed Hamiltonian 
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FIG. 1: Oscillator phase space diagrams for a system with a. = 2 and L = 4.74, and varying values of the energy i? = 0, 1, 2, 3, 6, 
and 10, as shown (See Sec. [H] for the definition of the reduced parameters used). The vertical axis in each figure is the scaled 
oscillator coordinate ( — d, and the horizontal axis the reduced oscillator momentum rj defined in Eq. 1181 



system. First, for any given system parameters, the phase space has a very simply described, completely integrable 
region of controllable size for all energies ranging from that of the ground state, Eg = —Fq/2Mlu'^, up to a critical 
positive energy Ea = \Eg\ (See Fig. n](a) and Fig.n(b)). For small positive energies the region of phase space outside 
of this region, which we refer to as Void I, appears in our numerical calculations to be fully chaotic, with no secondary 
KAM structures, even very close to the edge of the Void (See, e.g.. Fig. H^b)). This then makes for the simplest 
possible mixed phase space structure, in which motions in the completely integrable region and in the chaotic region 
are clearly separated and independent of one another. We are unaware of other Hamiltonian systems which exhibit 
this clean separation, although it does appear in a piecewise linear symplectic map on the torus that was explicitly 
constructed to exhibit this property 0. 

As a second striking feature of the model, we identify one-parameter families of marginally unstable periodic orbits, 
that appear as circular arcs in the oscillator phase space (See Fig. O, and are similar to the so-called bouncing ball 
orbits that arise in the stadium billiard. 

We furthermore show that for fixed, suitably-chosen system parameters, the statistical properties of the dynamics 
on the different energy surfaces can vary greatly. For some values of the energy, the dynamics is completely integrable, 
whereas for others it is fully chaotic or displays a mixed, KAM-type behavior. Among the rich variety of structures 
that arise, we identify and locate the central fixed point of a KAM-type elliptic island that we refer to as Void II. In 
some cases Void II appears alone, as in Fig.^d)-(f) and Fig. |2 while in others Void I and Void II both appear, as in 
Fig-EIb) andEIc). Finally, in the limit of small coupling strengths Fq, the phase space displays typical KAM structures 
of the type that generally occur when a completely integrable motion is subject to a small nonlinear perturbation, as 
in Fig. 121 This, however, by no means exhausts the variety of structures that seem to appear (See Fig. |SJ). 

The presence of chaos in our coupled particle-oscillator model can be understood intuitively as follows. A trajectory 
of the combined system will tend to be unstable whenever the time the particle takes between two impacts at x = ±(t 
is long compared to the oscillator period. Indeed, when the particle goes slowly, a small change in its velocity at one 
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FIG. 2: Oscillator phase space diagrams for a system with E = 0.85 and L = 2, and values of the coupling strength a — 0.7, 0.1, 
and 0.032 (See Sec.|n]for the definition of the reduced parameters used). The vertical axis in each figure is the scaled oscillator 
coordinate C, — d, and the horizontal axis the reduced oscillator momentum rj defined in Eq. 1131 The prominent elliptic island 
in these figures is what is referred to in the text as Void II. In the last two figures other secondary KAM structures appear at 
the edges of the Void. 

impact will lead to a large change in the relatively fast moving oscillator coordinate at the next impact. As a result, 
the height of the potential barrier the particle meets at that moment becomes highly unpredictable, and this is the 
source of the instability. That this simple picture is in agreement with observed behavior can be seen in Fig. |31 in 
which it is clear for a given _Fo and L (represented in that figure by the dimensionless parameters a and L introduced 
in the next section) , that the fraction of phase space outside Void I that is chaotic tends to increase as the total energy 
of the system (and thus the particle speed) decrease. Similarly, increasing L tends to increase the chaotic fraction of 
phase space outside of Void 1, since for larger L the particle spends more time between succesive visits at x = ±a 
whenever it is outside the interaction region. A more precise and quantitative analysis of the mechanism leading to 
chaos in this model will be given in Sec. IIIII where we will show that the dynamics can be analyzed in terms of a 
discontinuous linked twist map on the union of two intersecting disks. This is a generalization of the linked twist 
maps on the torus introduced and studied in Q. 

The rest of the paper is laid out as follows. In the next section, we simplify the Hamiltonian ^ by reducing the six 
system parameters with which it is associated down to two. Following this reduction, we present phase-space plots of 
the results of a numerical integration of the equations of motion for the system. In Sec. Illll we develop a geometric 
description of the dynamics that makes it relatively straightforward to explain the main features seen in numerical 
studies, including the emergence of chaos at positive energies (Sec. I1V|I . the existence of islands of regular motion of 
two different generic types, which we refer to as Void I and Void II (studied in Sees. IVl and I VII respectively), and 
the presence of arcs of marginally stable periodic orbits that appear as sets of zero measure in the chaotic portions 
of our phase space diagrams (Sec. 0). In addition, we analytically demonstrate the existence of unstable isolated 
periodic orbits of arbitrarily large period (Sec. IVl|l . In the last section, we summarize our results, and comment on 
their ramifications. 



II. DIMENSIONAL REDUCTION OF THE HAMILTONIAN 

The system described by the Hamiltonian _ff (x,po: P) as given in depends upon six system parameters: the 
two masses m and M, the oscillator frequency w, and the coupling strength Fq, as well as the widths 2a and L of the 
interacting and non-interacting sections of the ring on which the particle moves. To reduce the number of inessential 
parameters in the system, we now introduce a dimensionless time r = tot, and dimensionless variables and momenta 
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FIG. 3: Oscillator phase space diagrams for a system with a = 2 and L — 2.105, and energy E — 0.85, 1.35, and 3. The vertical 
axis in each figure is the scaled oscillator coordinate — d, and the horizontal axis the reduced oscillator momentum rj defined 
in Eq. [T31 



where h is an arbitrary constant having units of action that plays no part in the subsequent dynamics. The new 
variables obey the equations of motion 

dq . dp ^idx,.. 

^=n ^ = -$ + 5x(g) (3) 
which are the canonical equations generated by a transformed Hamiltonian — H/hw, where 

i/5 = ^(p' + n' + $')-5$x(9), (4) 
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a = ^ Fq /Alhuj-^, and the function x{q) = p(9\/fi-/wa;) vanishes outside the interaction region extending between 

1 /2 

q = ±(T = ±(T {rauj /h) . Explicitly, we write 

(5) 

where 6 {x) is the Heaviside step function. This form for the function x allows an additional simplification through 
the scale transformation 

q = q/(J p = 'p/(j $ = $/(T n = n/a. (6) 
With a suitable redefinition of the coupling constant a = aa, we obtain the following one parameter family 

Hc^\ {p" + n2 + $2) „ c.a>x(g) (7) 

of reduced Hamiltonians Ha = Ha<J~^ describing the system, where now the interaction region is associated with the 
fixed interval q £ [—1, 1]. Note that in this form the dimensionless coupling constant 



a^^^\E,\/]^ma^uJ^ (8) 

involves the ratio of the (original) ground state energy of the system to the kinetic energy of a particle that crosses 
the interaction region in one oscillator period. 

We have thus reduced the number of system parameters down to the single explicit parameter a describing the 
coupling strength, and one additional parameter L = L/a associated with the total range 2 -I- L of the particle 



1,0 
0.5 

-0.5 
-1,0 



1 1 1 1 1 


T = 2 - 


1 1 1 1 1 








'-■■v., (nl _ 
1 1 1 1 1 


1 1 1 1 


■■ (CI - 

1 1 1 1 1 



-1.0 



-0.5 0.5 1.0 -1.0 -0.5 0.5 1.0 -1 -0 5 0.0 5 1 



FIG. 4: Oscillator phase space diagrams for a system with E — 1.6 and a — 2, and values of the length L = 4,2, and 0.32, 
of the non-interacting region as shown (See Sec. [H] for the definition of the reduced parameters used). The vertical axis in 
each figure is the scaled oscillator coordinate ~ d, and the horizontal axis the reduced oscillator momentum rj defined in Eq. 
1131 At the center of each of these figures is a small Void I. As L is decreased, Void II appears and then ascends toward the 
impenetrable Void I. The collision between these two elliptic islands appears to generate considerable structure. 



coordinate q. We note in passing that with this choice for the function x (q) j the couphng parameter is equal to the 
equilibrium value <f>"q = a of the oscillator coordinate when the particle is in the interaction region q (z [— 1, 1] , and 
that the ground state, which now has rescaled energy Eg = — occurs when the particle and the oscillator are 
both at rest, with the particle in the interaction region. 
The equations of motion corresponding to ^ 

q=p p^a^[6{q + l)-S{q-l)] 

$ = n il= -^ + ax{q), (9) 

with dots denoting derivatives with respect to r, show that the particle feels an impulsive force only when it reaches 
the edges of the interaction region, but otherwise travels as a free particle. The particle's momentum undergoes 
discontinuous changes at these moments, but its position remains a continuous function of time. The impulse imparted 
to the particle at g = ±1 is readily computed from the oscillator displacement using only conservation of total energy. 
When the particle enters the interaction region, the oscillator experiences an interaction force of finite magnitude a 
that suddenly shifts the equilibrium position about which it oscillates from $°q' = to = a, but the amplitude and 
momentum of the oscillator remain continuous functions of time. Thus, when the particle is in the interaction region, 
the oscillator phase point (<I> (r) , 11 (r)) rotates at unit angular speed about the equilibrium state ($, IT) — {a, 0) , and 
when the particle is outside the interaction region, a similar free rotation takes place about the origin of the oscillator 
phase plane. As mentioned in the introduction, it is this essentially uncoupled evolution of the two subsystems 
between impulsive kicks of the particle at q = ±1 that makes it straightforward to numerically and analytically track 
the resulting dynamics. 

As in any system with two degrees of freedom, the energy surfaces associated with Q are three dimensional. Their 
two-dimensional sections at g = are ellipsoids centered at (<i>,n,p) = (a, 0,0). In Fig. |31we display phase points 
recorded on the Poincare section at g = and p > arising from the numerically determined evolution of 100 initial 
conditions randomly chosen on such an energy surface and evolved for a time corresponding to 500 passages through 
the section. Note that the Hamiltonian Q is invariant with respect to the parity operation of the particle. Thus, if 
(4> (t) ,U{t) ,p{t) ,q (t)) is a solution to the equations of motion, so is (<& (t) , 11 (t) , ~p (t) , —q {t)), and for symmetrized 
pairs of randomly chosen initial conditions, the corresponding set of phase points recorded at — p will be the same 
as at p, i.e., the back half of the ellipsoid, if displayed, would be a mirror image of the front half. This symmetry, 
together with time reversal invariance explains the additional 11 —11 symmetry that our phase plots exhibit. 

As a result it therefore suffices to represent the evolution by simply recording the oscillator phase point (<& (t) , tt (t)) 
each time that g = and p is positive. Such a representation is obtained by numerically computing the return map 
for this Poincare section, and has been used in Figs. 11141 and 1^ 
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FIG. 5: Visualization of the Poincare section at a; = of the energy surface for a particle-oscillator system with E = 1.6, 
Of = 2, and L = 0.32 (See Sec. ^ for the definition of the reduced parameters used). We display phase points generated, as 
described in the text, for 100 initial conditions randomly chosen on the energy surface. The same data appear in Fig. 2Jc). 
' ... - „ . , , , ated outside the enerffv surfa.ce. unside down in the TT-r) 




FIG. 6: The oscillator phase plane showing the disks D+{E) and D-{E) for (a) E < o? /2, and (b) for E > a^/2. On these 
figures the horizontal axis corresponds to the oscillator coordinate $ and the vertical axis to the oscillator momentum IT. 



III. GEOMETRIC DESCRIPTION OF THE DYNAMICS 

Although the dynamics of the system recorded at g = shows most clearly the inherent symmetries of the system, to 
actually explain the features that appear in our phase plots, it is obviously necessary to consider those times at which 
the particle reaches the edges of the interaction region at g = ±1. At fixed energy and with the position q = ±1 of 
the particle determined up to a sign, the state of the system at these moments is also conveniently represented as a 
point in the $ — 11 oscillator phase plane. In general, at a given E' > 0, the available phase space for the oscillator when 
the particle is in the interaction zone is a disk D^{E) of radius \/2E + cP- centered at the point (a, 0) in the ($,n) 
plane, that obviously contains the origin of the phase plane. Similarly, when the particle is outside the interaction 
zone, the available phase space for the oscillator is a disk D+{E) of radius \/2E centered at the origin. Any point 
(<I>,n) in the set Si-n (E) = D-{E) \ D+{E), i.e., inside D^{E) but outside £>+(£'), corresponds to a state with the 
particle definitely inside the interaction region, and any point in 5out {E) — -D+ {E) \ Z?_ [E) to a state with the 
particle definitely outside the interaction region. Depending on the energy and the coupling constant, two distinct 
situations can occur: either the disk {E) does not contain the center of D- (E), which occurs for E < 0^/2, and 
is depicted in Fig. Efa), or it does contain the center of D_ (E), which occurs when E > and is depicted in Fig. 
El^b). It should be obvious from the geometry of both figures that the circles which form the edges of the two disks 
always intersect on the 11 axis. In the following, given a point X = (<I>,n), we denote by A± its distance from the 
center of D±{E). Physically, = ^ + 11^) is the uncoupled oscillator energy. 
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We can now give a simple qualitative description of the motion of the system in geometric terms (see Fig. El^b)). 
Suppose that, for fixed E > 0, the system is such that at t = the particle is on the left edge of the interaction zone 
and moving inward, i.e., q{Q) = —1 and g(0) > 0, so that Xq — ($o,no) G D^{E). As the particle now crosses the 
interaction region to q — 1, the oscillator phase point simply rotates on a circular arc around the point (a, 0) through 
an angle that, in our dimensionless units, can be written 

-^(0) j2E + a^~A^ 



where the second form follows through energy conservation. Note that the rotation angle depends monotonically on 
the radius of the circle on which the phase point moves, so that two nearby phase points at different radii will rotate 
through different angles. In this way a shear T_ is induced on the disk (E) . Since the angle (p-{A^) diverges 
as A_ approaches the radius of D^{E), the strength of the shear diverges near the edge, and the dynamics becomes 
increasingly sensitive to small changes in A- . 

If the new phase point Xi so obtained lies in Sin (E) , the particle encounters a potential energy barrier at g = 1 
that is greater than its kinetic energy. Consequently, it reflects from the barrier, reversing its motion to cross the 
interaction region back to q = —1, during which time the oscillator phase point continues its rotation about the center 
of D-{E), through the same angle The process then repeats itself, generating a sequence of oscillator phase 

points Xi, X2 , . . . separated by equal angular displacements until, for some fc, the phase point Xk falls inside 

S{E). 

At such a time, the particle has sufficient kinetic energy to overcome the barrier it encounters, and passes out of the 
interaction region with a new kinetic energy equal to E — 5^^, where A+ is the distance between the point Xk and 
the origin (the center of the disk (E)). The particle then travels through a distance L around the outer section of 
the ring, and arrives again at q = ±1, while the oscillator phase point rotates on the second disk (E) through an 
angle 

^+(A+) = , ^ (11) 
2E-AI 

about the origin of the phase plane. Thus, a different shear T+ is induced on the disk £>+ (E). The new phase point 
Xk+i so obtained is depicted in Fig. |Sfb). Obviously this way of depicting the evolution can be repeated indefinitely, 
and provides an efficient way to analyze and explain various features of the dynamics. 



IV. THE EMERGENCE OF CHAOS 



For example, the description given in the last section makes it clear that the mechanism underlying the appearance 
of chaos in this system is the existence of the two non-aligned shears T_(_ and T_ . This is a well-known phenomenon. 
For example, on a torus x,y G [0, 1], successive application of the two shears 

yields a hyperbolic map 

{y)^{b l + ab){y) 

and leads to chaotic dynamics if |2 -I- ab\ > 2. Locally, the two shears r+ and T_ have exactly this structure, but in 
polar coordinates and with the role of a and b played by ip'^{A+) and (p'_{A-). The divergence of ip±{A±) at the 
edges of the disks thus provides a clear mechanism for the emergence of chaos in certain regions of phase space. In 
fact, the dynamical system defined on the two disks described in the last section is a generalization of what is referred 
to in the literature as a linked twist map. Some simple examples of such maps (on a torus rather than on a union of 
disks) have rigourously been shown to exhibit ergodicity and chaotic behavior [3 ■ The discontinuity of the functions 
(p± at the edges of the discs and the fact that the two shears are not transverse everywhere in the intersection region 
would make a completely rigorous analysis of the ergodic properties of our model considerably more complicated. 

Note that in this essentially geometric description, the dynamics depends only on three independent parameters 
E, L, and a. It is helpful to simplify the geometric nature of the description even further, by reorganizing the 



FIG. 7: The rcscalcd oscillator phase plane showing the disks D+ and The shaded region indicates those initial conditions 
ill D- which never enter the intersection S of D- and L>+, and thus load to trajectories in which the particle remains in the 
interaction region. These trajectories make up what we have referred to as Void I, with the center of D- corresponding to the 
center of that elliptic island. 



parameters as follows. First, we re-express the coupling strength of the system through the parameter 

d= " ^ : (12) 



and introduce rescaled oscillator variables 



$ n 

^lE^o?' V2E + a^ ^^^^ 



which locate the oscillator phase point at radii 



r± = -=4=, 0<r_<l, < r+ < - ^2 (14) 
V2E + 



from the centers, respectively, of a disk D- of unit radius centered at (d, 0) , and a disk of radius \/l — c?^ 

centered at the origin. With this choice, the dynamics now occurs on the disks D±, and the rotation angles (p± that 
the oscillator phase points sweep through during one traversal of the interaction zone can be simply written 

<p_(r_) = <^_(A_) = ^^, a. = ^=2== (15) 
/l _ J.2 ^2E + 



^+ {r+) = y+ {A+) = ^^1^ , a+ = ^===. (16) 



In this now description, wc can thus take d, a+, and a_ as the three independent parameters describing the system. 
Note that numerical data in the oscillator phase plots appearing throughout this paper are presented in terms of the 
scaled oscillator variables {( — d,rf) , and thus always appear on a disk of unit radius centered at the origin. This 
disk is essentially a shifted version of the disk D_ defined above, but is rotated by 90 degrees, so that the oscillator 
coordinate appears on the vertical, rather than the horizontal axis. In what follows we will use both the scaled and 
unsealed variables as is appropriate to the discussion at hand. 



V. MOTION CONFINED TO THE INTERACTION REGION 



In this section we analyze the motion of the system when the particle never leaves the interaction region. In 
particular, we explain the basic features of Void I, the region of complete integrability that occurs at the center of the 
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FIG. 8; A partial trajectory in the q — ^ plane for a motion in which the particle remains confined to the interaction region. 
For this system a — 1 and E — 0.2. Shaded portions indicate energetically inaccessible regions of configuration space for a 
system with this energy, and the region between the vertical dashed lines indicates the interaction region as described in the 
text. 

phase space plots in Figs.QJa) and (b), and in Fig. 01 In addition, we demonstrate and explain the existence of arcs 
of parabolic fixed points embedded in the chaotic regions of phase space. 

The easiest case of confined motion to understand is the one in which the total energy is negative, and the particle 
is energetically confined to the interaction region. For this situation the dynamics is totally integrable, since the speed 
and kinetic energy of the particle are constants of the motion, and the oscillator coordinate ^(t) is never negative. 
There are essentially three types of trajectories. The first type includes those in which only the motion of the particle 
is excited, so $ = a and 11 = are constant. Since the particle does not have enough kinetic energy to overcome the 
barriers it encounters at q = ±1, it reflects at each impact with the boundary. The trajectory is trivially periodic, 
with period A/q. In the second type of motion, the particle is at rest in the interaction region and the oscillator 
performs simple harmonic motion of amplitude < a about <i> = a. Periodic orbits of this type are also possible, 
of course, at any positive energy. Finally, there are confined motions in which both of the "modes" described above 
are excited. The trajectory is Lissajou-like, and the orbit is closed if the oscillator period and the traversal time 
2/q are commensurate. Otherwise it sweeps out a rectangle in the two-dimensional {q,^) configuration space, with 
q E [—1,1], and <i> g [a — ^o, a + $o]j for some positive amplitude <I>o < a of oscillation. 

When the total energy E — 0, possible motions with the particle confined to the interaction region include all of 
the types described above. In addition, oscillations of amplitude $0 = a are now allowed; in that case the particle 
is at rest somewhere in the interaction region. We note that that the only other motions at i5 = are the unstable 
equilibrium points corresponding to the particle at rest outside the interaction region with the oscillator at rest as 
well. 

We now give a description of all orbits at positive energy for which the particle never leaves the interaction zone. 
It is clear from Fig. |B|(a), the scaled version of which appears in Fig. [TJ and the geometric description of Sec. IIIII 
that for any energy > E > 0, the disk does not contain the center of Thus, for this situation, all initial 
conditions with r_ < i?i , indicated by the shaded region in Fig. [T] in which 

Ri = d- Vl - d^, (17) 

will give rise to trajectories for which the particle remains in the interaction region. Indeed, since the circle of radius 
r_ does not then intersect ZJ-i-, successive rotations of such a point through if - (r_) can never lead to a point lying 
in Z?+. The elliptic regions of complete integrability at the center of Figs. H^a) and (b), and in Fig. 01 which we 
have already referred to as Void I, correspond precisely to trajectories of this kind. In Fig. |H1 we show a typical 
periodic orbit of this type in the q — ^ plane. Shaded portions of that figure indicate classically forbidden regions 
at this energy. From the equivalent point of view of a particle moving in a two dimensional potential V {q, <&) , it is 
interesting that the motion remains trapped within the interaction region although there is no actual potential energy 
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FIG. 9: Phase plot showing arcs of periodic orbits. In this figure one can see a single arc of fixed points lying just outside 
of Void I, and, moving out from the center, sets of arcs corresponding to orbits of period fc = 3, 2, 3, and 2, corresponding, 
respectively, to ^ = 4, 3, 5, and 5. 



barrier preventing it from leaving. Finally, we note that as E is increased from zero at constant a the radius of Void I 
shrinks, as described by (|I7|I . until it finally disappears at the critical value E = (corresponding to c? = l/\/2), 

as in Fig.^c). For E <0 Void I fills the entire oscillator phase space, as in Fig.^a). 

Aside from these trajectories that remain within Void I, there are at any positive energy still other initial conditions 
that give rise to trajectories in which the particle remains in the interaction region. As we will show, however, such 
trajectories are then necessarily periodic. Indeed, suppose r_ > d — y^l — cP. Then the circle of radius r_ does 
necessarily intersect and successive rotations through (/3_(r_) will take the orbit into i?+ whenever the rotation 
angle (^_(r_) is an irrational multiple of 2n. Thus for the orbit to remain in Z?_, it is necessary that the rotation 
angle (^_(r_) = 2^1/ k be rational, and hence that the orbit be periodic. Note that k is then the period of the orbit 
(where it is understood that the integers k and £ are relatively prime). Furthermore, for such a fc— periodic orbit to 
be possible, the phase points must somehow arrange to miss the intersection with as they advance around £)_ in 
angular steps of 27r£/fc. Clearly, for this to happen, the angle 5 (r_) subtended by the two intersection points of the 
edge of the disk and the circle of radius r_ on which the phase point moves must be smaller than the angular 
separation between neighboring orbit points, i.e., 

5 (r_) < 27r/fc. 

For fixed points (fc = 1), this condition on the angle 5 (r_) is automatically satisfied, but for periodic orbits with 
fc > 2, the condition imposes a bound on r_ that prevents, e.g., periodic orbits of this type occuring too close to the 
edge of -D_. 

These arguments show that, in general, fixed points of the dynamics will occur at the radii 



f— V 



I > a_/27r 



(18) 



for which (fe,-) = 2111. Note that this infinite sequence of values r£^_ accumulates at r_ = 1, i.e., at the edge of 
D_. It follows that any phase point in — D^(E) \ D^{E) with r_ = r^__ will be a fixed point of the dynamics. 
The set of such points for a given value of I is an arc in the phase plane. In Fig. |31 such an arc of fixed points with 
I —1 appears near the edge of Void I. 

Now suppose for some value of fc > 1 and r_ we have (r_) = 27r£/fc, and 5 (r_) < 2ix jk. In this situation there 
will exist fc arcs of angular span 2iT/k — 5 (r_) centered at {d, 0) , each point of which is associated with a period fc 
orbit. Several arcs associated with higher order orbits of this type also appear in Fig.|^ and can occasionally be seen 
in some of our other phase plots. If the phase points were recorded at g = ±1, then each arc that appeared in such a 
diagram would, by the arguments given above, lie outside of the intersection region with Z?_, which is at the bottom 
of each phase diagram in this paper. However, because the phase points in these figures are recorded at g = when 
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p is positive, the angular position of each arc is rotated by an odd multiple of tp- (r_) /2 from where it would be if 
recorded at q = ±1. This is due to the rotation of the oscillator phase point that occurs while the particle travels 
from the edge of the interaction region back to the center, where the phase point is recorded. Since this rotation is a 
function of the radius r_, arcs of this kind can generally appear at any orientation in the phase diagram. 

On general arguments it is to be expected that periodic trajectories of this kind arc parabolic. Indeed, the stability 
matrix associated with a phase point lying on such an arc has a vanishing Lyapunov exponent along the direction 
tangent to the arc itself, since neighboring points along the arc are periodic orbits of the same order. On the other 
hand, in any neighborhood of such a point there will be initial conditions at radii r_ just above or below the arc that 
will not satisfy H18() for any values of t and k. Such points will give rise to trajectories in which the particle eventually 
does pass outside the interaction region, to end up at the mercy of the alternating shears and T_, which give rise 
to the chaotic portions of the phase diagram. Thus, in general, we expect arcs of periodic trajectories of this type to 
be largely immersed in the chaotic parts of the phase diagram. 

From the discussion above, which focuses on motions confined to the interaction region, one may wonder whether 
there are counterparts to these motions which take place with the particle confined entirely outside the interaction 
zone. In fact, it is easy to convince oneself that such motions are relatively few and far between. Indeed, for any 
point Xq in S'out, the circle on which it moves necessarily intersects Moreover, the angle (5+ (r_|_) subtended by 
the two intersection points of the edge of the disk Z?_ and the circle of radius r+ on which the phase point moves 
is now greater than tt. Hence, only period one orbits are able to avoid entering the interaction region. Such fixed 
points will occur whenever the time the particle takes to traverse the non-interaction region is an integer multiple of 
the oscillator period. 

VI. OTHER PERIODIC ORBITS: VOID II 

Having thus classified all orbits in which the particle never leaves the interaction region, we now turn to the more 
complicated situation in which the particle explores the entire configuration space available to it. Of course, it is 
impossible to classify all such orbits, since they are precisely the ones that are responsible for the variety of intricate 
structures that appear, as well as for the chaos. Nonetheless, some of the more prominent features of our figures can 
be explained quantitatively. Indeed, we show in the analysis below that the fixed point at the centers of the elliptic 
islands that we have referred to as Void II arise from a relatively simple type of orbit in which the particle traverses 
each section of the ring exactly once per period. We also show that an infinite number of such fixed points exists, and 
that most of them are hyperbolic, and thus unstable. 

Conceptually, such an orbit can be viewed within the geometric picture developed above, as follows. Consider a 
phase point Xq on the ^-axis within D_ at an instant when the particle is at the center of the interaction region and is 
moving to the right (See Fig. [rUT a)). When the particle reaches q = 1, the oscillator phase point has rotated through 
an angle ip- (r_) /2 about the center of D_ to the point Xi, as shown. The particle now leaves the interaction region 
and travels around the ring to q = —1, while the oscillator phase point rotates through an angle (r+) to a point 
X2. For the particle to now re-enter the interaction region, the point X2 must lie in the intersection region S, as 
shown in Fig. liur a'). If, in addition, the system is to return to its original state when the particle arrives again at 
g = 0, it is obviously necessary for X2 to fall on the circle of radius r_ on which it started, as in Fig. llOf b). Thus, 
all periodic orbits that traverse each section of the ring exactly once per period have vertices Xi and X2 that lie on 
a "lozenge" structure of this type. Using relatively simple arguments, given below, we can locate the tops of all such 
lozenge orbits and demonstrate in the process that they are infinite in number. 

To proceed, we recall that any point Xi which forms the top or bottom of a lozenge orbit can be be obtained by 
rotating some point on the segment (d — 1, d -I- 1) through an angle ip^ (r_) /2. The locus F^ of all such points is 
depicted in Fig. ^2 as a solid curve passing through the center of D- . But any such point Xi can also be obtained 
by rotating some point on the segment (— Vl — cP ^ Vl — <P^ backwards in time through an angle —p^ (r+) /2. The 
locus F+ of this set of points is depicted as a solid line passing through the center of D+ in the same figure. It should 
now be clear that the intersection points of these two curves locate all possible values Xi associated with periodic 
orbits of this type. For each Xi so located, the corresponding fixed point Xq is obtained by rotating Xi through 
—p>- /2. Since these two curves intersect an infinite number of times, the number of such fixed points is, itself, 

infinite. Indeed, the points and X^'^ at the intersections of the edges of D_ and D+ are accumulation points 

of lozenge tips. In addition, any point where either of the two curves plotted in Fig. ITTl crosses the edge of the disk 
in which it did not originate will also be an accumulation point of lozenge tips. As a result, it is clear that the points 
'^o±' with coordinates (1 ± d, 0) are accumulation points for fixed points of the dynamics. As we will show shortly, 
and as one might expect from the discussion of section HVI fixed points that occur too near the edge of either disk 
(where the shear strengths diverge) will be unstable. 

On the other hand, the fixed points at the centers of the elliptic islands visible in Figs, ^c)-(e), Figs. [2Ia) and 
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FIG. 10: Motion on the disks D- and D+, of lozenge orbits associated with fixed points of the dynamics. 



(c), Figs. |2Ic), and Fig. ^b) and (c) are clearly stable, and can therefore not be associated with an intersection of 
r+ and r_ that is too close to the edge of either disk. Based on this argument, the most reasonable candidate for 
such a fixed point is associated with the first intersection of r+ and r_ encountered when moving outward along r+ , 
starting from the origin. This intersection is readily computed numerically. For values of E^L, and a appropriate to 
our figures in which a Void II occurs, the results of such a calculation are tabulated in Tabled a-nd are indicated as 
horizontal dashed lines in the associated figures. To numerical accuracy the tabulated values agree with the actual 
locations of the centers of the type II Voids in those figures, supporting the basic picture developed above. 

Let us now show how to predict in advance whether any particular intersection of r+ and r_ gives rise to a stable or 
unstable fixed point. To this end, it is clearly sufficient to consider the stability of the evolution in the neighborhood 
of the upper lozenge tip. Let Xi and X2 be the vectors, respectively, from the origin to the upper tip Xi, and to the 
lower tip X2, of some lozenge, and let x'^ = xt — d be the corresponding vectors locating those points from the center 
of By definition, X2 = T^xi, and xi = T^X2, where T+ and T_ are the mappings associated with the evolution 
of the system during one traversal of the corresponding region by the particle. Note that T^T^xi — xi. It will 
therefore be sufficient to compute the Jacobian matrix of T-T^ at Xi. Because the points Xi and X2 are obtained 
from one another by moving along arcs of constant radius from the centers of D±, it follows that r± (xi) = r± {X2) , 
ip^ (xi) = iy9_(- {X2) = 0, and if - (xi) = ip^ (2^2) = V'- Specifically, this means that 

xi = d + R (tp) (^X2 — dj ~ r_ X2 X2 = R {((>) xi — Tj^xi 

where the rotation matrix 

f cose sine \ 
^ ' y—amO coso J 

induces a rotation in the clockwise sense through an angle 9. We compute first the Jacobian matrix J_|- of T+ at Xi 



J4 



cos(p sm(j) \ ^ (p / — sm0 cos0 \ / xiXi xiyi 
-sincf) cos(f> J r+ \-cos(j) -sint/) J y xiyi yiyi 

cos(j) sin0 \ (j)' f cos(0 — 7r/2) sin(0 — 7r/2)\ / xixi xiyi 
-sin0 COS0 J r+ \ - sin (0 - 7r/2) cos(0-7r/2) J \ xiyi yiyi 

= R{(t))~ r+(l)'R (0 - 7r/2) xiXi 

where xx denotes the second rank tensor product of the unit vector x with itself. Similarly, we consider the Jacobian 
J_ of r_ at X2 

J _ I cosip ainijj \ ^ ip' ( — sinV' cos-^ ^ ( x'^x'^ x'2y2 



-sinV' cosV' / r-{x2) V-costA -sinV' / \ X2y2 2/22/2 

cos^j) shxij} \ ill' f cos — 7r/2) sin (■(/; — 7r/2) \ f X2x'2 2:22/2 

-sinV' cosV' / r„ (^2) V - sin (?/; - 7r/2) cos (1/' - 7r/2) J \ X2y2 2/22/2 
R (tp) - r-ip'R {ip - 7r/2) £^^2. 
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FIG. 11: Intersection of the curves r_ and r+, described in the text, locating the tip Xi of a lozenge orbit. Each curve wraps 
an infinite number of times around the edge of the disk in which it originates. 



The trace of the product 

J_ J+ = (i? (V-) - r-i^'R {tjj - 7r/2) (0) - r+<j>'R {(/) - 7r/2) xixi) 

= R{4> + '4})- r-^'R {%p - 7r/2) x^x^R {(t>) - r+4>'R (V' + </) - 7r/2) xiXi 
+ r-r+ip'(t)'R (V- - 7r/2) x^x^R {(j) - n/2) xx 



of these two matrices determines the stabihty of the evolution in the neighborhood of the lozenge tip Xi associated 
with the orbit. The trace of the first term is obvious. To evaluate the second, use Tr [xy\ = x ■ y, and the invariance 
of the trace of a product of operators under their cyclic permutation: 

Tri? - 7r/2) x'^x'^R {(f) = TtR {(j)) R{ip - 7r/2) x'^x'^ = {R {(j) + - 7r/2) x'^) ■ x'^ 

= cos + — 7r/2) . 

Similarly, Ty [R {tp + (j) — 7r/2) xiXi] = cos {(j) + tp — tt/2) . The trace of the last term similarly reduces to the product 

Tr [R {iP - 7r/2) x'^x'^ R{(l>- 7r/2) x^Xi] = • i? (0 - tt/2) xi) Tr [R {%P - 7r/2) x'^xi] 

= (4 • [ii (.^ - 7r/2) xi]) (xi -[Rii^- 7r/2) x'^]) . 

To evaluate this we need to consider more carefully the geometry of the situation. Let 6- = ip/2 and 6+ = <f)/2 be 

the interior angles adjacent to the base of the triangle connecting the centers of _D_ and £)+ to the point Xi. Note 
that the sum 9+ + 9- < tt. The angle x = ti" — 6*+ — 0_ opposite this base is the angle between the vectors xi and 
xi — d = x'l- So 

xi=R{tt-9+ -9^)x[. 

But a vector along £'2 is rotated into a vector along x[ by a rotation through (j) = 29^, i.e., x[ = R (26*-) ij, so 



Xi 

This leads to 

xi-[R{'il)-n/2)x'2]=x-i- R 

and 

x'2 ■ [R{(j)- tt/2) xi] = x'2 ■ 
= cos 



R{ti - 9+ - 9-) R{29-) x'2 =R{tt-9+ 



.)x'2. 



f + ^-3./2)f, 



= cos 



= — sm 



i^ + cj) 



R{<j)-TT/2)R{T: -9+ + 



R 



+ I + 7r/2 ) x'2 



2 + i+^/2 



■ sm 
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Figure Co - d Tr[J-J+] 
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TABLE I: Location and stability of some lozenge type fixed points Xq — (Co,??o) = (Co,0) appearing in the figures indicated. 



Thus, for the lozenge evolution, 



Tr [J_J+] = 2cos(i/' + 0) - (r_V' + f-+'/'')sin(0 + V') +r-r+V'Vsm2 ( "^^-^ ] (19) 



Tr [J-J+] = 2 + (r_r+Vy0' - 4) sin^ ( ^^-—] - 2 (r_ V + r+0') sin ( ^^^-^ ) cos ( ^^-^ ] . (20) 



From this form it is first of all possible to prove, confirming our previous arguments, that fixed points associated with 
lozenge tips close to the edge of either disk (where ip' and diverge) are unstable. This is easily shown for d < 1/^/2. 
In that case, a little trigonometry shows that the angle x associated with any lozenge tip near the boundary of S is 
acute, so that tt > {ip + (j>) /2 > n/2. Consequently, Tr [J- J+] > 2, with equality only if sin {ip + (/)) /2 — 0, so that 
the fixed point is hyperbolic or, at worst, parabolic. For this last case to occur the angles -0/2 and 0/2 must add up 
to TT, which can only happen if r+ and r_ intersect on the horizontal axis. In that case, (j>/2 = n and -0/2 = 0, or vice 
versa. It is not hard to convince oneself that, for any fixed value of d, there exist arbitrarily large values of a_ and 
a+ where such fixed points will appear. All other fixed point of this type with d < 1/V2 are, however, hyperbolic. 

Using the algorithm described above we have calculated numerically the location and stability, as indicated by 
the value of Tr [J_ J_(.], for the first intersection of r_|_ and r_ for the system parameters of a certain number of our 
figures. The results appear in Table HJ Comparison with the relevant figures shows that whenever |Tr [J_J+] | < 2 
the associated fixed point is indeed the center of a Void 11. For Fig. Hb), this fixed point, computed to lie at 
— d — —0.575, has a trace with magnitude just greater than 2. It appears as the hyperbolic structure located below 
Void 1 in that figure. In Fig.^b), we have, in addition, computed the second and third intersection of r+ and r_. 
One leads to a fixed point at ( — d = 0.878, which is stable, and appears at the center of a crescent shaped strucure 
at the top of that figure. The other occurs at — d = —0.864, but is unstable. It does not, therefore, give rise to an 
elliptic island, but lies right at the edge of Void 11. 

Finally, although in general the location of the fixed point at the center of Void 11 must be determined numerically, 
its behavior for large energy can be obtained from the limiting behavior of the dynamics as d, a+, and a_ go to zero. 
In this limit, denoting the location of the fixed point as Xq = {rjd, 0) , we note that the curves T± are well represented 
by straight lines in the neighborhood of the origin. The corresponding triangle bounded by the segment [0, d] and 
r± then has a height h — rjdtana-i- ~ iida+ — rjLd^/a as measured from the origin, and h = {1 — rj) dtana- ~ 
(1 — 77) da- = (1 — 77) d^/a as measured from {d, 0) , so that 



locates the corresponding fixed point on the segment [0, d] . This has a simple physical interpretation. At high energies, 
when the particle is moving very fast, the fractional change in its speed as it enters and exits the interaction region is 
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small. In such a motion, the reduced oscillator coordinate oscillates about equilibrium position d during that fraction 
2/ (2 + L) of the time the particle is in the interaction region, and oscillates about the origin during that fraction 
L/ {2 + L) of the time that the particle is outside the interaction region. A time average of these two values of the 
equilibrium position results in the location H21I) of the fixed point in this high energy limit. 

Having thus made our point that the most prominent features that appear in our phase plots can be explained 
in terms of the simplest periodic orbits of the system, we point out that many other periodic orbits occur in which 
the particle traverses the two sections of the ring more than once per period. Such orbits will then give rise to 
the more complicated structures appearing, e.g. in Fig. [S] We conjecture that a more complete analysis of those 
orbits along lines similar to those developed in this paper, would allow a complete explanation of some of the more 
anthropomorphic structures appearing in that figure. 

VII. DISCUSSION AND SUMMARY 

The model we have introduced and studied in this paper clearly has certain similarities to a number of previously 
studied dynamical systems, such as the kicked rotor, the Fermi accelerator, billiards, and the spring-pendulum. Indeed, 
the study of the dynamics reduces in all these cases to the study of a return map on a suitably chosen Poincare section. 
The model presented here nevertheless differs in important ways from each of these previously studied systems. 

Indeed, the kicked rotor and the Fermi accelerator describe externally perturbed nonconservative systems with one 
degree of freedom, while the current model is closed, energy-conserving, and has two degrees of freedom. In our model, 
although the particle gets "kicked" each time it reaches q — ±1, the kicks are neither periodic in time, as in the rotor, 
nor arc they imposed by an external agent, as is the case for both the rotor and the Fermi accelerator. 

Billiards, of course, are closed conservative systems with two degrees of freedom as well, but trajectories in billiards 
have the unusual property of being independent of particle energy, so that changing the energy does not change the 
statistical features of the dynamics. This is in sharp contrast to the behavior of the present model, in which many of 
the interesting features that arise, do so as a result of changes in the energy of the system under conditions in which 
the underlying potential is kept fixed. Moreover, the present system can be interpreted in terms of the Hamiltonian 
interaction between two otherwise separate mechanical systems. It should, therefore, be of more use in understanding 
many problems for which that is an essential feature. 

Another closed Hamiltonian system of that type is the spring-pendulum Like the current model, the spring- 

pendulum is a Hamiltonian system with two degrees of freedom, each of which has a simple description when treated 
on its own. Unlike the spring-pendulum, however, the present model has the advantage that the coupling between the 
two sub-systems can be smoothly turned off in a way that allows the unperturbed dynamics of each to be recovered, 
and thus may be more useful for understanding fundamental properties of interacting independent systems of this 
type. We note that a bifurcation analysis for the spring-pendulum was given in j^. A similar analysis applied to our 
model would yield predictions on the motion and the shape of the elliptic islands that appear in our model as the 
parameters are changed, and could be an interesting direction for future study. Our emphasis here has instead been 
on some unusual features of our model, that we now briefly recall. 

We have argued that for suitable system parameters the system exhibits either a fully chaotic phase space, or a 
mixed phase space in which only two regions occur. In one the motion is chaotic, and in the other it is completely 
integrable with no secondary KAM structures (See Fig.^b) and (c)). While we have explained the presence of chaos 
in the model in terms of alternating shears similar to the kind that arise in linked twist maps, our conjecture regarding 
the absence of secondary structures of finite measure in the chaotic sea for small energies is based largely on numerical 
calculations over very small regions of phase space that have heretofore failed to detect any structure in the chaotic 
region. We have not, however, given a rigorous proof that such structures do not emerge at small enough length scales 
in phase space. 

We note also that the sharp boundary that exists in this model between the chaotic and the completely integrable 
parts of phase space should make it an interesting system on which to test current conjectures of quantum chaos 
theory, in particular those pertaining to systems with a mixed phase space, and to the localization properties of 
eigenf unctions on chaotic and completely i nteg rable parts of phase space. While the clear-cut boundary of the present 
model should facilitate such an analysis, 6, 10] it does have the complication of an underlying potential that has finite 
discontinuities in configuration space, and the matching problems that occur, in contrast, e.g., to billiard systems, 
for which a number of efhcient numerical techniques have been developed for solving the corresponding Schrodinger 
equation ITU . 

Our own interest in the present system arose originally from a fundamental interest in Hamiltonian models of 
transport and dissipation in deformable media. In many such systems, the medium through which a particular 
transport species moves can be modeled as an appropriate collection of harmonic oscillators. Among the many 
questions that arise in such extended systems is, e.g., to what extent the presence of microscopic chaos, or its absence, 



16 



in the interaction of a particle with a single oscillator, manifests itself in the transport properties that emerge at 
long times after repeated interactions with many independent or mutually- coupled oscillators. We view the present 
analysis as a step towards answering this and other questions related to systems of this kind. 
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